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Abstract. Germanium dioxide (GeC>2) is a chemical analogue of SiOi. Furthermore, 
it is also to some extent a structural analogue, as the low and high-pressure short-range 
order (tetrahedral and octahedral) is the same. However, a number of differences exist. 
For example, the GeOi phase diagram exhibits a smaller number of polymorphs, and all 
three GeOi phases (crystalline, glass, liquid) have an increased sensitivity to pressure, 
undergoing pressure induced changes at much lower pressures than their equivalent 
SiOi analogues. In addition, differences exist in GeOi glass in the medium range 
order, resulting in the glass transition temperature of germania being much lower 
than for silica. This review highlights the structure of amorphous GeOi by different 
experimental (e.g., Raman and NMR spectroscopy, neutron and x-ray diffraction) and 
theoretical methods (e.g., classical molecular dynamics, ab initio calculations). It also 
addresses the structure of liquid and crystalline GeOi that have received much less 
attention. Furthermore, we compare and contrast the structural differences between 
GeOi and SiOi, as well as, along the GeOi ~ SiOi join. It is probably a very timely 
review as interest in this compound, that can be investigated in the liquid state at 
relatively low temperatures and pressures, continues to increase. 
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1. Introduction 



Zachariasen P proposed the continuous Random network model (CRN) to explain the 
structure of oxide glasses, and it has subsequently received wide acceptance in describing 
glasses that form continuous random networks. To date the majority of studies of oxide 
glasses have involved the investigation of silica (SiOz) or borate (B2O3) glasses with, to 
a lesser extent, germania {GeO^) glasses. The structure of the latter has generally been 
considered to be comparable to that of silica glass despite differences in bond lengths, 
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angles and the relative size of Ge versus Si. Experimental studies of amorphous Ge0 2 
have generally involved either x-ray or neutron scattering and spectroscopic techniques 
such as x-ray absorption spectroscopy (EXAFS/XANES) and Raman spectroscopy. 
On the other hand, theoretical studies have generally employed classical or ab initio 
molecular dynamics calculations to gain insight into the structure of these materials. In 
both approaches, the results of the studies are often compared to the known crystalline 
polymorphs of Ge0 2 . Here we review the structure of amorphous Ge0 2 (glass and 
liquid) from both an experimental and theoretical perspective, as well as, comparing 
their structure with that of amorphous Si0 2 (glass, liquid). Furthermore, we review the 
structure of the crystalline Ge0 2 polymorphs, both at room temperature and pressure 
and at elevated temperatures and pressures. 

2. Crystalline Ge0 2 Polymorphs 
2.1. Structure 

Crystalline Ge0 2 exists at ambient temperatures and pressures as one of two polymorphs 
(Figure 1): an a-quartz-like (P3 2 21) trigonal (hexagonal) structure [2] or as a rutile-like 
tetragonal (P4 2 /mnm) structure 

The a-quartz-like Ge0 2 structure has been shown to be the stable high temperature 
phase jl] and, while the structure is very similar to that of a-quartz, there are some 
distinct differences. In particular the GeO^ tetrahedra are more distorted due to greater 
variation in the O — Ge — O angles within the tetrahedron, which range from 106.3° 
to 113.1° with a Ge — O — Ge angle of 130.1°. This is in contrast to a-quartz where 
the O — Si — O angles within the S2O4 tetrahedron are relatively uniform ranging 
from 108.3° to 110.7° with a Si - O - Si angle of 144.0° 5J. These differences 
are important as it results in different mechanisms being responsible for the high 
pressure behaviour of a-quartz and a-quartz-like Ge0 2 . For a-quartz the tetrahedra 
are relatively rigid and compression of the structure occurs via cooperative rotation or 
tilting of the tetrahedra around the shared bridging oxygens. Conversely, for a-quartz- 
like Ge0 2 while compression does occur via tilting of the tetrahedra, distortion of each 
tetrahedron via changes in the individual O — Ge — O angles also plays a large role 

This behaviour is the reason why germanate analogues of silicate phases are useful 
in high pressure studies since they undergo pressure induced phase transformations at 
much lower pressures than their silicate analogues. The transformation of a-quartz to 
the high pressure rutile structure (Stishovite) occurs at 10 GPa while the equivalent 
transformation for a-quartz-like Ge0 2 to rutile-like Ge0 2 has been observed to occur 
at much lower pressures when the sample is heated: ~ 1.8-2.2 GPa at 417 K 

As noted above, the stable room temperature Ge0 2 phase is the rutile Ge0 2 
polymorph which transforms to the a-quartz-like structure at 1281 K ([3], [2j and 
references therein). The rutile Ge0 2 polymorph has a structure similar to that of 
stishovite [3] and, like stishovite, the 2 axial bonds within the GeO§ polyhedron are 
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Figure 1. Projection of the a-quartz-like structure (left) and rutilc-like structure 
(right) onto the (001) plane. 

longer than the 4 equatorial Ge - O bonds: 1.902 ± 0.001 A and 1.872 ± 0.001 A, 
respectively. Conversely, the two independent Ge — O distances in the a-quartz-like 
Ge0 2 structure are similar at 1.737 ± 0.003 A and 1.741 ± 0.002 A 0. 

2.2. High pressure and temperature behaviour 

A number of studies have investigated the high pressure behaviour of the two Ge0 2 
polymorphs. Itie et al. [Hj investigated a-quartz-like Ge0 2 at ambient temperature. 
They observed an increase in the Ge — O bond length and Ge coordination number 
consistent with the formation of the rutile-like Ge0 2 phase between 7—9 GPa. However, 
subsequent studies have suggested that the transformation is to an amorphous phase 
rather than the crystalline rutile-like Ge0 2 polymorph jjj, 0, JU|, CB- Furthermore, 
it has been suggested that the amorphization step is a precursor to subsequent 
transformation to the rutile polymorph jTTj. More recently, Brazhkin et al. El 
have shown that with compression, a-Ge0 2 changes via a martensitic transition into 
a crystalline monoclinic (P2i/c) phase. On the other hand, Haines et al. P3J suggest 
that there is no evidence for amorphization of the crystal. Instead a poorly crystalline 
monoclinic (P2i/c) phase forms consisting of edge sharing chains of GeO^ octahedra 
(Figure 2). 

The monoclinic phase is metastable up to 50 GPa. However, when combined with 
heating, it transforms to the rutile structure at pressures up to 22 GPa and above 
43 GPa forms a mixture of CsC^-type and i^iV-type (or a-Pb0 2 , see later) high 
pressure phases [TJ] , jTHj , [15] • This monoclinic phase was also reported by Prakapenka et 
al. JHl between 7 — 52 GPa at room temperature but with laser heating it transforms to 
an orthorhombic CaC/2-type structure above 36.4 GPa and an a — Pb0 2 -tjpe structure 
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Figure 2. Polyhedral representation of the (3 x 2)-kinked P2\/c structure of GeC>2 
determined by Haines et al. |15|, Crystallographic axes are not to scale and are merely 
to show orientation of the structure. Green (light) octahedra are the fully occupied Ge 
positions while the blue (dark) octahedra are the Ge octahedra which exhibit partial 
occupancy. 

at 41 GPa. 

Haines et al. have also observed transformation of the rutile-like GeO^ phase 
to the orthorhombic CaC7 2 -type structure above 25 GPa at ambient temperature while 
Ono et al. ^H] observed the transition at high pressure and temperature. Ah initio 
calculations by Lodziana et al. [TH] suggested that rutile-type GeO<i should transform 
to a-P&CVtype (above ~ 36 GPa) and pyrite (Pa) type (above ~ 65.5 GPa) structures 
and these were subsequently observed by Ono et al. [20], [21] around 44 and 90 GPa, 
respectively. A Pe2iV-type (or defect Ni— As) phase at pressures larger than 25 GPa has 
been observed by Liu et al. [22] and Haines et al. [T^j. This type of structure is similar to 
an a-P&CVtype structure but with the Ge sites disordered and has, more recently, been 
explicitly identified by Ono et al. [2H] as being the a-P&CVtype structure. However, 
it should be noted that Prakapenka et al. [2B] observe the defect NiAs structure when 
amorphous GeOi is heated to 1000 — 1300 K at 6 GPa (see later). 

Structural refinements of the crystalline phases have been obtained by Shiraki et al. 
[21] and a phase diagram for crystalline Ge02 is given in Figure 3. In addition, another 
orthorhombic phase has been suggested to occur at ~ 28 GPa and 1273 K by Ming and 
Manghnani [23]. They concluded that this phase was not the a-P&CVtype structure 
but it has not been observed subsequently. The phase transformation sequence of rutile- 
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Figure 3. Phase diagram of crystalline GeOi (after [23]). 

like Ge0 2 (25 GPa) CaCl 2 -type -> (44 GPa) -> a-P60 2 -type -> (70 - 90 GPa) 
pyrite-type structures is consistent with the high-pressure behaviour of other group-IV 
element dioxides such as Pb0 2 , Sn0 2 , and Si0 2 (cf., [26], [TH], [TO] , |23J ) . 

Of interest is the way in which the a-quartz-like and rutile-like Ge0 2 structures 
respond to increasing pressure. As noted above, Jorgensen |5J observed that in the a- 
quartz-like polymorph compression up to 2.5 GPa occurs predominantly via changes in 
the individual O — Ge — O angles and that tilting of tetrahedra was secondary. Yamanaka 
and Ogata [27] carried out a series of structural refinements on the a-quartz-like Ge0 2 
polymorph up to 4.48 GPa and found that the GeO^ tetrahedra are relatively rigid 
with little change in the Ge — O bond length, consistent with the study of Itie et al. 
[Hj. Yamanaka and Ogata .27J found that the dominant mechanism responsible for the 
observed pressure induced unit-cell volume change in the structure involved a decrease 
in the Ge — O — Ge angle from 130° to 125°. Conversely, Glinneman et al. j2H] found that 
tetrahedral tilting was responsible for the 11% volume change of a-quartz-like Ge0 2 up 
to 5.57 GPa. 

The phase transformation of rutile-like Ge0 2 to the CaCl 2 -type structure occurs via 
compression of the axial Ge — O bonds of the octahedron. The axial bonds are elongated 
relative to the equatorial bonds (see above). With increasing pressure there is increased 
compression of the axial relative to equatorial bonds and the transformation at 
25 GPa occurs during flattening of the octahedra [21]. With transformation to 
the a-Pb0 2 -tjpe structure, the GeO& octahedron becomes further deformed with the 
Ge atom displaced from the center of the octahedron and 2 of the six Ge — O bonds 
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becoming elongated [23], as suggested by the numerical results of Lodziana et al. |19J. 
Transformation to the pyrite-type structure, however, results in GeO§ octahedra that 
are symmetrical with Ge in the centre. 

As noted above, the a-quartz-like polymorph is the stable high temperature phase 
and rutile-type Ge02 will transform to this polymorph above 1320 K\ the transformation 
temperature being the highest of any of the quartz-like analogues. The high temperature 
(up to 1344 K) behaviour of this polymorph has been investigated by Haines et al. "21"] 
who found that the intertetrahedral bridging angle {Ge — O — Ge) and tilt angles exhibit 
thermal stabilities that are amongst the highest observed for quartz-type analogues. 
With increasing temperature, expansion of the unit cell is highly anisotropic with 
expansion along a being 5 times greater than along c "2""]. However, the a-quartz-like 
Ge02 polymorph is metastable at low temperatures ([30 provide a number of methods 
for growing the a-quartz-like Ge02 polymorph) but does undergo transformation to 
the rutile-type polymorph at around 1000 K although the reaction proceeds slowly due 
to the kinetics involved (cf. [HI])- Finally, it should be noted that a cristobalite-like 
polymorph for GeO<i has been observed after long-time heating of Ge02 glass to 873 K 
[3*2"] or by dehydration of ammonium hydrogen germanate ((NHf)3HGejOiQAH20) 
between 853 — 873 K [53], however, this polymorph has not been observed in in-situ high 
pressure and temperature studies. In addition, the a-quartz-like polymorph at 1322 K 
mentioned by Leadbetter and Wright [""3" and Desa et al. t35. based on the work of 
Laubengayer and Morton [3] and Sarver and Hummel "3""] has also not been observed. 

3. Ge02 glass structure 

3.1. Neutron and X-ray diffraction 

Neutron and X-ray diffraction data are complementary tools for inferring structural 
information since the chemical sensitivity is different for the two techniques; Ge — O 
and Ge — Ge pairs are better resolved with X-ray and Ge — O and O — O with neutron. 

Ge0 2 glass structure has been studied using X-ray diffraction in the pioneering 
work of Warren [HZ], [HE] and Zarzycki [21]], [HI]- It was found that the Ge atoms are 
arranged in basic tetrahedral units such as those found in the trigonal a-quartz-like Ge02 
polymorph. X-ray diffraction data with higher real space resolution (Q m ax = 17 A -1 ) 
confirmed these findings [33] and determined the first Ge — O and Ge — Ge distances 
at 1.74 A and 3.18 A, respectively, giving an intertetrahedral angle of ~ 133°. The 
first neutron diffraction experiment (Q m ax = 18 A^ 1 ) on vitreous Ge02 shows two 
strong peaks at 1.72 A and 2.85 A ascribed to Ge — O and O — O correlations, which 
is consistent with GeO^ tetrahedra [33] • The Ge — Ge peak, initially determined at 
3.45 A [31], [32], was resolved in a high resolution neutron diffraction investigation 
{Qmax = 35.5 A -1 ) at 3.21 A, which is slightly higher than the Ge — Ge distance 
determined by X-ray diffraction due to the overlapping of Ge — O and — pairs 
[4*5j . [33] • A recent neutron and X-ray diffraction investigation "3j"] has shown that the 
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O — Ge — O intratetrahedral angle is more distorted in vitreous Ge0 2 than in vitreous 
Si0 2 , with a distribution likely comparable to that of Ge0 2 a-quartz (106.3 — 113.1°). 
This is due to the larger radius of Ge than Si, allowing more accessible positions for O 
atoms around Ge atoms. The mean Ge-O-Ge intertetrahedral angle was estimated from 
the Ge-0 and Ge-Ge distances to be 130.1° with a range of 121 — 147°. This mean value 
was confirmed at 133 ± 8.3° using high-energy X-ray diffraction 03] • This bond angle 
and its distribution are lower than in the case of vitreous silica. The smaller Ge-O-Ge 
angle probably results from the presence of increased numbers of 3-membered rings in 
the Ge0 2 network relative to vitreous Si0 2 (cf. [US] and see later) since such planar 
rings have a Ge-O-Ge angle of 130.5° 0E],[E]. The values for the main interatomic 
distances, coordination numbers and intertetrahedral angles found in these studies are 
reported in Table I. The structure of Ge0 2 can thus be viewed as continuous random 
network of corner sharing tetrahedral, as in silica, but with greater distortion of the 
tetrahedra and larger amounts of three-membered rings. 

The diffraction data (Fig. 4) of Ge0 2 are composed of three partial functions, Ge- 
Ge, Ge-0 and O-O. The first attempt to separate the three components was carried out 
using X-ray anomalous diffraction and neutron diffraction 0H|, |49j . The Ge-O, 0-0 
and Ge-Ge distances are found at 1.73, 2.85, and 3.17 ± 0.04 A, respectively, and the 
average Ge-O-Ge intertetrahedral bond angle is estimated to have values between 129° 
and 139°. Recently, by combining neutron and X-ray diffraction, together with X-ray 
anomalous scattering, the three partial functions were fully separated up to Q = 9 A -1 
[ED! j EH] j though problems exist due to different instrumental resolution functions that 
appear especially at low Q values, and the necessity to improve the anomalous scattering 
terms [H2j. The structure factors are dominated by peaks occurring at 1.54, 2.6 and 
~ 4.5A" 1 . The first feature at 1.54 A -1 (usually called first sharp diffraction peak, 
FSDP) corresponds to intermediate range ordering and is stronger in X-ray than in 
neutron data (SBj- The FSDP is associated with a positive peak in SceGe and Sc e o 
and a shallow negative peak in S001 indicating that cation correlations dominate the 
medium range order (Figure 5, jSI])- The peak at 2.6 A -1 is strong and positive in 
SceGe and SoO) and strong and negative in Sceo, and has been attributed to chemical 
short range order. The peak at 4.4 A -1 occurs predominantly in SceGe and is due to 
topological short range order. The latter conclusion is usually extracted from Bhatia- 
Thornton structure factors [53| that show the correlations between number density and 
concentration fluctuations [HE]- On this basis, the chemical and topological ordering in 
Ge0 2 can be rationalized in terms of an interplay between the relative importance of 
two lengthscales that exist in the glass [56] . 

There have been considerable efforts to compare diffraction data obtained on Ge0 2 
glass with equivalent calculations based on the Ge0 2 crystalline polymorphs, with 
divergent results. Leadbetter and Wright [SI] concluded that the intermediate range 
order in the glass closely resembles a quasi- crystalline model based on the a-quartz-like 
Ge0 2 structure with a correlation length of 10.5 A but discrepancies appear beyond 
4 A. Bondot jinj obtained good agreement between the glass and the a- and /3-quartz 




Figure 4. Measured total structure factor (circles, [57]) together with calculated S(Q) 
from ab initio (solid line, |58|) and Classical Molecular Dynamics (broken line, |59|). 



Ge02 polymorphs, which led to the conclusion that the glass contains six-membered 
rings. On the contrary, Konnert et al., j^I] concluded that vitreous germania, like 
vitreous silica, possesses the same short range order as that found in the tridymite SiC>2 
polymorph. The vitreous GeC>2 structure could thus be described as randomly oriented, 
slightly distorted tridymite-like regions having dimensions ranging up to at least 20A 
\<nl\ . However, these regions are not crystallographically ordered (i.e. not microcystals) 
but have similar bonding topology in the glass and in tridymite. In a more recent 
study it was shown that, though similarities exists with crystalline a-quartz and 
a-cristobalite Ge02 polymorphs, diffraction data are not consistent with large volume 
fractions of quasi crystalline- like regions, due to an important distribution of torsion 
angles. 

3.2. Neutron and X-ray diffraction at high pressure and temperature 

Due to the technical difficulties associated with performing in situ diffraction 
experiments, pressure effects have been mainly studied on pressure-released glasses, 
in which permanent densification is observed. Permanently densified glasses (up to 
18 GPa) were studied by X-ray diffraction in the low Q-region (FSDP) which is sensitive 
to medium range order A shift to higher Q and an increase in width of the FSDP 
is observed above 6 GPa, a pressure corresponding to the threshold for coordination 
changes observed in in situ experiments (see below). However, comments on this study 
pointed out that changes in the diffraction peaks may not necessarily be associated with 
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Table 1. Interatomic distances (R), Coordination numbers (N), standard deviations 
(a) and Ge-O-Ge intertetrahedral angle determined by diffraction methods. 1) 
Qmax = 50 A" 1 , EZI; 2) Q ma:c = 50 A"\ 3) Q max = 9 A" 1 , EU; 4) 
[SSI; 5) 1 A" 1 < Q < 10 A" 1 , |5D1; 6) 0.6 A" 1 < Q < 33.5 A" 1 , gS]; 7) 
0.22 A" 1 < Q < 23.6 A" 1 , 8) [H|; 9) 0.8 A" 1 < Q < 17 A" 1 , ND 
= Neutron diffraction; AXS = anomalous X-ray scattering; HEXRD = high energy 
X-ray ray diffraction; XRD = X-ray diffraction. * Fixed values. 
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Q [A 1 ] * (A) 

Figure 5. Left from top to bottom: Measured (circles, [ST]) and calculated (solid 
line, [SOI) partial structure factors, SceGeiQ), Sceo(Q) and Soo(Q)- Right : Partial 
correlation functions, Qmm(t), gMo{ r ) and goo( r ) (from top to bottom), for the three 
atom pairs in vitreous GeOi at room temperature (M = Ge, points), together with 
the corresponding functions from rescaled Molecular Dynamics simulation of vitreous 
Si0 2 (M = Si, lines) (After [ET]). 

a coordination change [S3]. 

Ge0 2 glasses densified up to 6 GPa at 673 K (densification of 16%) were 
investigated by neutron diffraction [2Z|, while a glass densified at 10 GPa and 300 K 
(densification of 11%) was studied by neutron and X-ray diffraction (Figure 6) |55] . No 
evidence of six-coordinated Ge was observed. The GeO^ tetrahedra are distorted, with 
Ge — O distances increasing by 0.005 ± 0.001 A and O — O and Ge — Ge distances 
decreasing by 0.023 ± 0.002 A and 0.019 ± 0.002 A, respectively jH^J. The main change 
is a shift of the Ge — Ge peak (at ~ 3.1 A) to lower r values with increasing pressure 
compaction [HZj- This indicates a reduction in the mean Ge — O — Ge bond angle with 
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increasing density. Noticeable changes are seen for the FSDP in the neutron and X-ray 
structure factors: the FSDP shifts towards higher Q, broadens and become less intense 
on densification. This indicates a reduction of the network connectivity. By combining 
neutron and X-ray diffraction up to Q = 30 A -1 , it was shown j^Sl that variations of the 
FSDP are mostly associated with O — O correlations rather than Ge — Ge ones. This is 
attributed to a decrease in the average size of the network cages (these can be considered 
as holes in the structures, formed for instance, by the ring structures), yielding better 
packing of the GeO^ tetrahedra. 

In situ measurements have recently been obtained (Figure 6) by both neutron (up 
to 5 GPa) and X-ray (up to 15 GPa) diffraction [66 . The FSDP decreases and almost 
vanishes with increasing pressure in neutron measurements while it gradually shifts 
to higher Q in X-ray data. This is interpreted as a breakdown of the intermediate 
range order upon compaction of the tetrahedral network associated with changes in 
the oxygen correlations. In the X-ray correlation functions, a reduction of the Ge — O 
distance is observed below 6 GPa while it increases at further pressure, corresponding 
to GeO^ tetrahedra being converted to GeO^ octahedra. Based on molecular dynamics 
simulations, it was argued that stable five-fold units are present in the transition 
region, indicating a new intermediate form of the glass. The structure of the high 
pressure glass is based on edge- and corner-shared octahedra, which is not retained 
upon decompression. 

Ge0 2 in the liquid state has been investigated by X-ray diffraction [10], jHZj. The 
Ge — O distance is unchanged in agreement with a small thermal expansion of the Ge — O 
bond similar to that for Si — O bonds. The GeO^ tetrahedra are preserved in the Ge0 2 
melt but Ge — Ge distances are shifted from 3.16 A at room temperature to 3.25 A at 
1100°c, which is interpreted as a widening of the Ge — O — Ge bond angle. 

3.3. Raman Spectroscopy 

3.3.1. Ge0 2 polymorphs The Raman spectra of the crystalline polymorphs of Ge0 2 
(Figure 7) were first reported by Scott |68J. The rutile-like Ge0 2 spectrum exhibits 
three strong bands in the 150 — 1200 cm -1 range at 173, 701 and 873 cm" 1 . The band 
at 701 cm" 1 is the Ai g mode while the 873 cm" 1 band is the B 2g mode. The Bi g mode 
is at 173 cm" 1 . The E g mode observed at 680 cm" 1 by Scott [HE] is not observed in the 
spectrum shown in figure 7. Alpha-quartz-like Ge0 2 has a number of additional bands 
including four symmetric modes of Ai symmetry and 8 doubly degenerate modes of E 
symmetry all split into transverse optic (TO) and longitudinal optic modes (LO) [68J. 
The a-quartz-like Ge0 2 spectrum of figure 7 is comparable to that first obtained by 
Scott jSHJ. Bands can be assigned following Scott [68J and Dultz et al. [HH] as A 1 modes 
at 263, 330, 444, and 881 cm" 1 ; E modes at 123 (TO+LO), 166 (TO+LO), 212 (TO), 
330 (TO), 516 (LO), 593 (LO), 860 (TO), 960 (TO), and 973 cm" 1 (LO). E modes 
at 372 (LO), 385 (TO), 492 (TO), 583 (TO) and 949 cm" 1 (LO) are too weak to be 
observed in the spectrum or are unlabelled for clarity. 
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Figure 6. In situ structure factors (after jHEj for (A) X-ray diffraction at 0, 3, 5, 6, 
7, 10, and 15 GPa (bottom to top) and (B) neutron diffraction up to 5 GPa, with 
ambient-pressure data from Sampath et al. |f>5| . 
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Figure 7. Unpolarizcd Raman spectra of the trigonal and rutile GeOi polymorphs. 
Bands are comparable to those of Scott [HE] and Madon et al. [HI]. For clarity not all 
bands are labelled. 
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The effects of increasing pressure and temperature on the vibrational spectra of 
the Ge0 2 polymorphs has been investigated by Sharma [70], Madon et al. [3T] and 
Mernagh and Liu [H] • With increasing pressure the Raman bands for the a-quartz-like 
polymorph shift to higher wavenumber and behave in a similar manner as the IR modes 
[31]. The mean frequency shift is ~ 1 cmT^/O.l GPa for bands in 400 — 600 cm~ l region, 
0.3 cm^ 1 /O.I GPa for the bands in the 100 — 330 cm -1 region and the bands in the 
850 — 970 cmT 1 region do not shift at all up to 4 GPa except for the band at 961 arc 1 . 
The rutile Raman bands behave slightly differently jZO] with the band at ~ 173 cmT x 
shifting to lower wavenumbers and the other two bands to higher wavenumbers. With 
increasing temperature, the Raman bands of the rutile-like polymorph transform to the 
a-quartz-like spectrum at ~ 1313 K [31] while the Raman bands of the a-quartz-like 
polymorph show a nonlinear shift with increasing T. Madon et al. [Hi] observed a shift 
of —0.01 cm -1 / K for the bands in the low frequency region and —0.024 cm~ 1 /K in the 
mid- and high-frequency regions. For the rutile polymorph, the Raman bands above 
600 cm~ l exhibit nonlinear shifts to lower wavenumbers whereas the 173 cm -1 band 
exhibits a shift to higher wavenumbers with increasing T |71j . In addition Mernagh and 
Liu [7_T| detect (by deconvolution) splitting of the A\ g mode (701 cm -1 ) of the rutile-like 
polymorph with a new band observed at 684 cm -1 . 

3.3.2. Ge0 2 glass and liquid The first Raman spectrum of Ge0 2 glass was described 
by Bobovich and Tolub [Z2] and Obikhov-Denisov et al. [23]. A Raman spectrum for 
Ge0 2 glass is shown in figure 8a. The Raman band assignments for Ge0 2 glass are 
similar to those of Si0 2 glass but are shifted to lower frequencies (wavenumbers, cm~ l ) 
because of the larger mass of Ge relative to Si. Currently accepted band assignments for 
Ge0 2 are given in Table 2 and extensive discussion of Raman assignments and earlier 
literature can be found in [23], |7K] . 

The high frequency bands observed at ~ 860 and 998 arc 1 are the TO and LO split 
asymmetric stretching bands of the bridging oxygens (Ge — O — Ge). The Ge — O — Ge 
bending modes are observed in the broad region between ~ 500 — 620 cm -1 and have 
also been assigned to TO (~ 556 cm -1 ) and LO (~ 595 cm -1 ) split modes associated 
with significant Ge and O motion [ZOj. A " defecf band D 2 occurs at ~ 520 cm -1 . This 
defect mode is the equivalent of the D 2 band observed in Si0 2 glass at 606 cm -1 and 
is assigned similarly to an oxygen-breathing mode associated with 3-membered rings 
of Ge0 4 tetrahedra and this assignment has been recently supported by the study of 
Giacomazzi et al. [58J (Fig. 8c). The intensity of this band is much stronger relative 
to the main vibrational band at ~ 420 cm -1 in comparison to the equivalent bands for 
Si0 2 glass. This indicates that the Ge0 2 network, while being composed predominantly 
of 6-membered rings of GeO^ tetrahedra (see above) does have a larger proportion of 3- 
membered rings relative to Si0 2 glass. The relatively narrow band at around 420 cm -1 
is the symmetric stretching mode of the Ge — O — Ge bridging oxygens. Its width is much 
narrower than the equivalent band observed in Si0 2 glass at 440 cm~ l and indicates 
that the distribution of Ge — O — Ge intertetrahedral angles for Ge0 2 glass is narrower 
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Figure 8. a) Unpolarised Raman spectrum of GeOi glass showing the main vibrational 
bands, b) a Raman spectrum of a L^O-containing germanate glass showing the high 
frequency BO and NBO bands: The insert is a curve fit (deconvolution) of the high 
frequency envelope into its discrete vibrational bands (see table 2), c) Unpolarised 
Raman spectrum of GeOi glass (solid line) compared with calculated HH spectrum of 
Ge0 2 glass (dashed-dot line) from Giacomazzi et al. The calculated spectrum 

has been shifted so that the main vibrational band is coincident with the equivalent 
band of the experimental spectrum. 
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Frequency 


Attribution 


~ 60 cvrT x 


Boson peak; Acoustic mode? Related to glass fragility 


o4< cm (.-UIJ 


Ge " deformation" motion within the network 


420 cm" 1 


Symmetric stretching of bridging oxygens (BO) in 6-membered 
GeOi rings: Ge-O-Ge 


520 cm- 1 (D2) 


" Defecf mode assigned to breathing motion of bridging 
oxygens m o-memueieu ^orew^ imgs 


OUU — DZU C777. 


Bending modes: Ge-O-Ge, TO (556 cm ) and LO split 
(595 cm' 1 ) 


~ 860 cmT 1 


TO asymmetric stretching of bridging oxygens: Ge-O-Ge (Q 4 ) 


r^j O Q Q ™^ — 1 

— yoo C777, 


LO asymmetric stretching of bridging oxygens: Ge-O-Ge (Q 4 ) 


Q species vibrations 


NBO (non-bridging oxygen) vibrations that occur upon 
addition of network modifiers such as alkalis 


~ 865 cmr 1 


Q 3 Ge0 4 tetrahedra with 3 BO and 1 NBO 


780 cm- 1 


Q 2 Ge0 4 tetrahedra with 2 BO and 2 NBO 



Table 2. Raman band assignments for Ge02 glass and for the different Q species 
observed upon addition of a network modifier such as an alkali or alkaline-earth cation. 



than that for S2O2 glass consistent with the neutron and X-ray data above. 

The origin of the Boson peak (BP) at 60 cm' 1 (the peak occurs over a broad range 
between 40 — 60 cm- 1 ) remains controversial. It has been assigned to acoustic-like 
harmonic modes, localized quasi-harmonic modes, and to the smallest energy van Hove 
singularity of the crystal (cf. [77], [7H] and references therein). Most recently there 
seems to be a consensus that the origin of the BP is due to optic-like excitations related 
to nearly rigid SiO^ (or GeO±) librations through hybridization of the acoustic waves 
[TU] . It exhibits a dependence on the fragility of the glass (fragile glasses have weak BP 
intensity), as well as, fictive temperature (for S1O2 the BP shifts to higher wavenumber 
with increasing fictive temperature). In addition, there is a monotonic frequency shift 
in the BP for pure S1O2, to lower wavenumbers with the addition of Ge02 which may 
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indicate that Ge0 2 substitutes isomorphously into Si0 2 [7H] . 

The Raman spectrum of Ge0 2 glass indicates that the T-O-T intertetrahedral 
angle and its distribution are narrower for Ge0 2 glass relative to Si0 2 glass, consistent 
with the X-ray and neutron diffraction studies (see above). These latter studies also 
suggest that the medium-range structure of Ge0 2 glass consists of 6-membered rings 
of Ge04 tetrahedra, similar to those observed in the a-quartz polymorph of Ge0 2 , 
with a high proportion of small 3-membered GeO± rings (relative to Si0 2 glass). An 
interesting aspect of the medium-range structure was raised by Henderson et al. [13] and 
Henderson and Fleet [75] using Raman spectroscopy. They suggested that the medium- 
range structure of Ge0 2 may actually consist of 4- rather than 6-membered GeO^ rings. 
This suggestion has not been explored further and there have not been any X-ray or 
neutron scattering studies that have compared Ge0 2 glass with structures containing 
predominantly 4-membered GeO^ rings. However, Giacomazzi et al. [HE] recently used 
a model Ge0 2 structure that had exclusively 3 and 4-membered GeO± rings. Their 
model reproduced the first sharp diffraction peak in the neutron static structure factor 
(indicative of medium range structure), and the infrared and Raman spectra of Ge0 2 
glass (Figure 8c) reasonably well. The question of whether or not the medium-range 
structure of Ge0 2 glass consists of 6- or 4-membered rings remains unanswered and 
open for further studies. 

3.4- Infra Red (IR) spectroscopy 

There have been relatively few infra-red studies of Ge0 2 glass primarily because the IR 
spectra are more difficult to interpret and obtain than the Raman spectra. One of the 
earliest is that of Kaiser et al. jHO] while more recent studies have tended to use IR in 
high-pressure studies jSIj for investigating the onset of amorphous to amorphous phase 
transitions (see below). The IR spectrum of Ge0 2 glass exhibits two peaks at 560 crrT x 
and one at ~ 870 cm~ l with a shoulder at ~ 1000 cm -1 although the relative intensities 
for these two bands are reversed in the spectra of Galeener et al. [52] ■ The low frequency 
band at 560 cm -1 is the IR equivalent of the LO bending mode observed in the Raman 
spectrum at ~ 595 cmT x while the bands at 870 cm -1 and ~ 1000 cmT x are the IR 
equivalent TO (870 cm -1 ) and LO split asymmetric stretching of the bridging oxygens 
|82j . The data of Galeener et al. [82J also show a peak in their IR reflectance spectrum 
at ~ 340 cm -1 which is the equivalent of the 347 cm -1 Raman band. Galeener et al. 
[8*2] assign this band however to an LO mode. In general the LO modes are more intense 
in the IR relative to Raman spectra while the TO modes are more intense in the Raman 
relative to IR spectra. 

With increasing pressure, the 560 and 870 cm -1 peaks broaden and the region 
between the bands (~ 700 cm -1 ) exhibits an increase in intensity [SI], although part 
of this increase is due to a shift in the 560 cm -1 band to higher wavenumbers with 
increasing pressure (up to 6 GPa). Teredesai et al. [HI] also observe with increasing 
pressure a decrease in wavenumber for both high wavenumber bands. Above 6 GPa, 
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all bands shift to higher wavenumbers coincident with the onset of the pressure induced 
coordination change of Ge noted by Itie et al. jHJ. However glasses decompressed from 
9.5 GPa exhibit a 30 cm -1 red shift in the position of the 870 cm -1 peak with no shift 
in position of the 560 cm -1 peak |HT] . 

3.5. Increasing Pressure and Temperature 

The effect of pressure on GeC>2 glass at ambient temperature has been investigated by 
Ishihara et al. 83; and in situ by Durben and Wolf |84 and Polsky et al. [85.;. Up to 
6 GPa, Durben and Wolf [HI] observe a shift of the main Raman band at ~ 420 cm" 1 
to higher frequency with concomitant broadening and loss of intensity. Between 6 and 
13 GPa the main Raman band broadens and losses intensity without a shift in frequency. 
In addition they observe the growth of a broad low frequency band at ~ 240 cm -1 and 
no further spectral changes are observed beyond 13 GPa up to 56 GPa. However, upon 
decompression, the 520 cm" 1 D2 band characteristic of 3-membered rings is enhanced 
relative to uncompressed Ge02 glass and indicates that 3-membered rings are formed 
during decompression from high pressure. Similar results were obtained by Polsky et 
al. JH2] and both they and Durben and Wolf [84J observe subtle changes in the Raman 
spectra between 5 and 10 GPa characteristic of the pressure induced change in Ge 
coordination observed by in-situ EXAFS and XANES studies jSj. 

Up to 5 GPa, both Durban and Wolf [El] and Polsky et al. [ES] suggest that 
compression of the GeO^ glass network is taken up be tetrahedral deformation with a 
smaller decrease in the intertetrahedral angle. In addition, they conclude that there is no 
increase in intensity of the 520 cm~ l D2 band. However, this conclusion is questionable 
given that the main Raman band at 420 cm -1 appears to move to higher wavenumbers 
with increasing pressure (above 4 GPa) and as Polsky et al. [S3] themselves note, any 
apparent decrease in the intensity of the 520 cm. -1 band may simply be a consequence 
of changes in the adjacent band at 420 cm" 1 . Examination of Figure 2 of Durben 
and Wolf [84J shows that above 3.7 GPa the 420 and 520 cm -1 bands are merged and 
individual bands are unable to be discriminated. Below 3.7 GPa, the intensity of the D 2 
band also cannot be determined without some knowledge of how the spectra have been 
normalized but a cursory examination appears to indicate that the D2 intensity has 
increased relative to the maximum in the main 420 cm -1 band. Furthermore, Ishihara 
et al. [HH], albeit using permanently densified Ge02 glasses, note that growth of the 
D 2 band correlates with increasing pressure; higher pressures produce increased D 2 
intensity although there are no permanent structural changes for glasses decompressed 
from below 4 GPa 

High temperature studies have been performed by Magruder III et al. [SH] and 
Sharma et al. [HZ]- With increasing temperature Magruder III et al. determined that 
the high frequency TO/LO split pair undergo a 2-fold loss of intensity between 1723 and 
2023 K and that the D2 band intensity remains constant. The high frequency LO band 
at 988 cm -1 (figure 7a) loses intensity as the TO/LO splitting is lost with increasing 



CONTENTS 



19 



temperature [S7j but even in the melt phase two bands are observed at ~ 818 and 
940 cm -1 , respectively. However, Sharma et al. jHZl observed an increase in intensity of 
the D2 band and a shift of the main Raman band at 420 cm -1 to higher wavenumbers 
combined with a loss of intensity while the low frequency band at 347 cm -1 shifts to 
lower wavenumbers but is observed up to 1623 K. Both studies clearly show that the 
Raman bands observed in GeO^ glass remain even in to the melt phase but that there 
are subtle changes in intensities and band positions as the glass is heated and eventually 
melts. 

3.6. NMR spectroscopy 

The coordination environment of Ge in Ge02 and alkali-containing Ge02 glasses remains 
an area of intense interest from a glass perspective because of the unusual physical 
properties of alkali-containing germanate glasses and the possible role of Ge coordination 
in this behaviour [ZE], |88j . Ge NMR would normally be the technique of choice to 
investigate the coordination environment of Ge in glasses. Germanium has five naturally 
occurring isotopes ( 70 Ge, 72 Ge, 73 Ge, 74 Ge and 76 Ge) but only 73 Ge is suitable for NMR 
studies. However, while 73 Ge NMR has been successfully performed on solid crystalline 
compounds jEHj, |90] . [9T] it has not been useful for elucidating the structure of glasses 

The 17 O MAS NMR spectra of GeO<i glass, and the a-quartz-like and rutile-like 
polymorphs of crystalline Ge02 have been obtained by Du and Stebbins (OBI- The 
two crystalline polymorphs and Ge02 glass all exhibit a single crystallographic oxygen 
site similar to previous data obtained at lower magnetic fields [01]. The oxygen site 
in the Ge02 glass is comparable to that found in the a-quartz-like Ge02 polymorph 
indicating that the glass consists of a network of GeO^ tetrahedra, consistent with X- 
ray and neutron scattering studies. 

4. Structure of densified liquid Ge0 2 

Melting curves at elevated pressures were first reported by Jackson in the range 
1100 — 1700°C and 0.5 — 2 GPa. The high-temperature part of the phase diagram was 
also studied in [96J where an observed flatening of the melting curve at P ~ 2 — 4 GPa 
seems to be an indication of densification of the melt due to the transformation of a 
quartz-like liquid into a rutile-like one. 

Ordering of the melt structure in the same range of temperatures and pressures as 
above was also reported |HZj from Raman scattering. Specifically, the lowering of the 
Raiyleigh line intensity from in situ high pressure and temperature liquid Raman spectra 
was found to be significantly lower than for a glass quenched at ambient pressure. This 
suggest an increased degree of short-range order on compression in the liquid and a more 
ordered network structure. However, it constitutes a major obstacle to studying liquid 
Ge02 at elevated pressures. Note that Moelcular Dynamics has not tried to simulate 
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these experiments yet (see below). 

5. Structure of the binary Si0 2 -Ge0 2 glasses 

Germania and silica are prototype glasses for continuous random network models, based 
on the corner sharing connection of their SiO± and GeO^ tetrahedra. The variations 
in the intertetrahedral angles and the presence of some structural defects (for instance 
dangling bonds in SiC>2 glass) allows the formation of a three dimensional disordered 
network. 

Germanosilicate glasses are widely used as low-attenuation optical fibers, yielding 
numerous studies on their physical (optical) properties [OS]- Structural studies are 
more scarce despite the need for an understanding of the relationship between glass 
properties and structure, particularly with respect to variation of the local site geometry, 
intertetrahedral angles, ring statistics and their relationship to chemical ordering, 
clustering and/or substitution. A fundamental question is to determine whether or 
not germanosilicate glasses form a homogeneous network or if there is some sort of 
clustering or phase separation. 

5.1. EXAFS and X-ray scattering 

An early Ge K-edge EXAFS investigation 99 on 12.5Ge0 2 -87.5Si0 2 and 36.5Ge0 2 - 
63.5Si0 2 glasses calculated a Ge — O distance of 1.73 ± 0.01 A but no second 
neighbors were observed. A more extensive study using a combination of Ge K-edge 
X-ray absorption and wide angle X-ray scattering (WAXS) experiments were carried 
out on Ge0 2 — Si0 2 glasses containing 16 — 36 mol% Ge0 2 |l()()j . They showed that 
the XANES spectra are similar with increasing Ge0 2 content and that EXAFS-derived 
distances are 1.72 ± 0.02 A for Ge — O. A Ge coordination number of 3.9 ± 0.2, 
consistent with Ge in tetrahedral sites as in vitreous Ge0 2 , and a mean Si—O distance of 
1.62 A consistent with Si remaining tetrahedrally coordinated, were obtained from their 
WAXS data. These results seem to be contradicted by a high energy X-ray diffraction 
study on a 29Ge0 2 — 71Si0 2 composition glass [10 1] that found a mean coordination 
number for Si and Ge of 3.4 ± 0.05. The authors explain this low coordination number 
by proposing that a considerable number of Ge atoms are connected with less than four 
oxygens or are highly distorted. Except for the latter study whose coordination number 
seems questionable, all structural studies are consistent with the presence of SiO^ and 
Ge04 tetrahedra in binary Si0 2 — Ge0 2 glasses. 

The first peak observed in the X-ray radial distribution function is at higher distance 
than would be expected assuming standard Si — O and Ge — O distance (1.62 and 1.72 A 
respectively for tetrahedral environment |100j ). This suggests that the binary glasses 
are not a simple physical mixture of Si0 2 and Ge0 2 oxides. In germanosilicate glasses, 
no Ge0 2 clusters are observed and Ge04 tetrahedra are thus part of the Si0 2 network. 
This is confirmed by the second shell of neighbors that has been observed in EXAFS data 
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|lUUj . Indeed, this peak corresponds to Si and/or Ge neighbors and both its position 
and its intensity vary upon Si/Ge substitution. Ge atoms can thus be accommodated 
within the Si0 2 network. Intertetrahedral angles were calculated from EXAFS and 
WAXS data and are between 139 — 149°, which is closer to 144° for Si0 2 glass (albeit 
the magnitude of this angle remains controversial, cf. |l()2j than 133° for Ge02 glass. 
This suggests that at low Ge02 content, the Ge environment is constrained by the 
silicate network. These results are consistent with a substitutional model in which Ge 
substitutes randomly for Si in the vitreous Si02 network with little Ge clustering. A 
random substitution model is further supported by recent 17 O multiple quantum NMR 
spectra on Ge0 2 — Si02 binary glasses which show peaks for all three types of bridging 
oxygens (Ge — O — Ge, Ge — O — Si, Si — O — Si), in proportions at least roughly 
consistent with random mixing of the tetrahedral cations |103| . 

The binary Si0 2 — Ge0 2 glass structure can be described by a continuous random 
network of corner sharing GeO^ and SiO^ tetrahedra. 

5.2. Raman spectroscopy 

Information at the medium range structure such as ring statistics and the ordering of 
Si and Ge atoms, have been primarily obtained by Raman spectroscopic investigations 
of germanosilicate glasses. Important modifications appear between the Raman spectra 
of pure Ge0 2 and Si0 2 and some specific structures are present in the spectra of the 
binary glasses. 

The band at low frequency shifts from 437 cm~ l in Si0 2 to 416 cm -1 in Ge0 2 and 
becomes sharper |l()4j . This band is attributed to the T-O-T (T=Si or Ge) symmetric 
stretching mode and is thus characteristic of the distribution maximum in the T-O-T 
intertetrahedral angles |lU5j . Therefore, it can be concluded that the fluctuation in the 
intertetrahedral angle decreases as Ge0 2 is introduced into the silica network. 

A complete Raman study from pure Si0 2 to pure Ge0 2 was carried out by Sharma 
et al. jl()4j in order to characterize the distribution of SiO± and GeO^ tetrahedra. 
In the germanosilicate glasses, a weak band in the range 970 — 1010 cmT x appears 
that is not present in pure Si0 2 or Ge0 2 glasses. This band is attributed to the 
antisymmetric stretching motion of the bridging oxygen of Si — O — Ge linkages, while 
the corresponding modes for the Si — O — Si and Ge — O — Ge linkages appear at 
~ 1110 cm -1 and ~ 880 cm -1 , respectively. The position of the band is at ~ 1000 cm -1 
for the 10GeO 2 - 90SiO 2 glass but decreases to ~ 920 cm" 1 for the 90GeO 2 - !0SiO 2 
glass. This shift in position towards lower frequency is attributed to a decrease in the 
Si — O — Ge bond angle in the Ge0 2 -rich glasses |104j . In the 50GeO 2 — 50SiO 2 glass, 
the bands at 1100 and 880 cm" 1 are stronger than the one at 980 cm -1 . This indicates 
the formation of Si — O — Ge bonds but also the existence of an important number of 
Si — O — Si and Ge — O — Ge linkages. According to these authors, the Si/Ge ordering 
is likely non-ideal, which supports a random distribution of SiO^ and GeO^ tetrahedra 
(see comment above regarding 17 O NMR). In a Molecular Dynamics simulation of a 
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50GeO 2 — 505*202 glass |106| . a large fraction of Ge — O — Si bonds were found, as well 
as, Ge — O — Ge and Si — O — Si linkages. Based on the simulations, Bernard et al. 106 
proposed that Ge/ Si ordering occurred but not to the extent that phase separation was 
evident. They also showed that non-bridging oxygens (5%) were mainly localized in the 
Ge environment. 

With a small addition of Ge02, the Di and D 2 lines of Si0 2 glass at 495 
and 606 cm" 1 , attributed to four-membered and three- membered rings of the SiO^ 
tetrahedra in vitreous silica, are still observed but the intensity of the D\ line decreases 
sharply, while that of the D 2 lines decreases slowly and broadens |107j . Nian et al. |107j 
suggested that the substitution of Ge for Si in the vitreous Si02 network prevents the 
formation of these ring structures. This was explained by the disruption of the fourfold 
and threefold SiO^ rings to accommodate the larger GeO\ tetrahedra that distort the 
silicate network. Alternatively, the decrease in intensity of D\ and D 2 lines could also be 
due to a change in polarizability of the Si — O bonds as Ge pulls electron density away 
from O attached to Si, which is an explanation more consistent with the preference for 
3-membered rings in Ge0 2 relative to Si0 2 as indicated by diffraction and Raman data 
(see above). With small addition of Ge0 2 , a new band appears at 710 cm" 1 but its 
assignment is not clear |107j . Above 15 mol% Ge0 2 content, weak and broad shoulders 
are exhibited at ~ 568 and ~ 670 cm" 1 and the band at ~ 800 cm," 1 decreases in 
intensity |104j . 

5. 3. Evolution with pressure 

The Ge coordination change in the tetrahedral framework Si0 2 — Ge0 2 glasses is a 
reversible process that has to been studied by in situ high-pressure XAS measurements 
at the Ge K-edge |108j . The pressure-composition diagram in Fig. 9 shows the existence 
of three regions with distinct short-range structures. At low pressure, the region 
corresponds to a tetrahedral framework structure (T domain), then, an intermediate 
domain with a mixture of different sites, while, at higher pressure, the Oc region 
corresponds to a structure with ^Ge. The ^Ge to ^Ge transformation is a reversible 
with an important hysteresis (a return back to the tetrahedral site below 4 GPa). The 
coordination change is dependent on the mean composition of the glasses and extends 
over higher-pressure range when the Si0 2 content increases. The remarkable dependence 
of the Ge coordination change on the Si0 2 content shows that the Ge local structure 
is strongly affected by Si. The disruption of the Si0 2 tetrahedral network begins at 
10 GPa, as evaluated by Raman spectroscopy which is similar to the XAS data at 
high content and indicates that Ge and Si convert to a sixfold coordination state 
simultaneously. This result suggests that the pressured-induced transformations occur 
homogeneously in the mixed network and may be driven by the oxygen atoms rather 
than by the Ge or Si atoms [8*5] . Indeed, in such fully polymerized networks, oxygens 
increase their coordination from two to three in the transformation. The transformation 
occurs at higher pressure and over a broader pressure range when the Si0 2 content 
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Figure 9. Pressure-composition diagram of the <S«02 — GeOi glasses (after |1()8| ) 
depicting the pressure-induced Ge coordination change. The horizontal dashed line 
separates the pressure-range explored. Diamonds are the pressure onsets of the 
coordination change and the dashed curves delimit the intermediate domain separating 
the low-density form (T for tetrahedral) and high-density form (Oc for octahedral). 



increases. A careful analysis of both XANES and EXAFS signals supports a model of a 
mixing of ^Ge and ^Ge states in the transition region, in agreement with a kinetically 
hindered first-order process for the transformation at room temperature. 

6. Molecular simulations and theoretical approaches 

6.1. Force field parameters 

Several force field potentials have been proposed to describe either the crystalline phases 
of GeOi or amorphous germania. All these potentials contain a long range Coulombic 
part, along with a short-range repulsive term and an additional van der Waals-like term 

V ij (r ij ) = ^ + A ij e-^-^ (1) 

where Z± is the charge on ion %. The parameter serve to determine the steepness 
of the short range repulsive potential and is known as the "softness" parameter. The 
parameters and Cy serve to adjust the positions of the first peak in each possible 
radial distribution function to experimental findings. Oeffher and Elliot [109] have fitted 
equ. (UJ) to obtain cell parameters, density and elastic constants of the trigonal a-quartz 
like and tetragonal rutile-like phases of Ge02- Bond angles and bond lengths in both the 
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low and high pressure phases are found to agree with experimental findings. The Raman 
and Infrared vibrational spectra are also simulated within the harmonic approximation 
using the bond-polarizability model of Long Analysis from the vibrational density 

of states of Ge — O — Ge motions shows that for a-quartz-like GeO^-, symmetric and 
asymmetric bending motions are mostly confined to medium and low frequency bands 
while symmetric stretching and bending motions can be reasonably simulated at the 
anticipated frequencies. 

Matsui and co-workers have used the same kind of approach, i.e. the fitting of 
equ. ((TJ), to simulate another structural phase transition, namely the pressure induced 
change from a-quartz-like Ge02 to rutile-like Ge02 which happens at 7.4 GPa. The 
structure obtained at this pressure appears to be quite similar to the structure calculated 
for S1O2 at 21.5 GPa |112j . Furthermore it is shown that a-quartz-like Ge02 close 
to the transition is mechanically unstable as some of the elastic moduli of the lattice 
become negative. Specifically, the decrease of the transverse elastic constant C44 leads 
to an unstable shear that originates the transformation to the rutile-like structure. 
For increased pressures, a post-rutile-like structure is found ^Jj that has a CaC/2-like 
structure which consists of tilted GeO^ octahedra. This appears to be in agreement 
with Brillouin scattering results of a — Ge02 under pressure |113j which show that the 
shear constants are largely softened with respect to Si0 2 and can be related to shear 
instability. 

More recently, an alternative model has been proposed by Van Hoang |114j for 
liquid and amorphous germania that is based on a Morse-like potential in a similar 
manner to the potential given by Kim for Ge0 2 |115| . We discuss below the structural 
predictions of the Van Hoang potential. For completeness, we mention also the model 
potential proposed by Nanba |116j to account for Ge0 2 — PbO — PbF 2 glasses. However, 
it appears to show poor agreement with the rutile-like properties of GeC^- 

Topological and geometrical approaches have also been proposed |117j in order to 
generate continuous random network models of Ge02 that reproduce the experimental 
density, bond angle distributions and neutron scattering data [35J. Araujo has used 
statistical mechanical techniques |f 18j to calculate the density of oxygen vacancies in 
Ge02 and the absorption coefficient with respect to temperature. 

6.2. Simulation of liquid and amorphous germania 

Most of the work using the effective potentials described above, has been devoted to the 
description of the high temperature liquid where experimental data is lacking. Gutierrez 
and Rogan |119j have simulated Ge0 2 at 1500 K and 3000 K. At these temperatures, 
the system seems to be made of slightly distorted GeO^ tetrahedra which are linked 
by corners and have a Ge — O — Ge angle of 130°, similar to the experimental value 
in the amorphous phase (Ge02 glass). A volume collapse, in the range 4 — 8 GPa, is 
seen from the pressure-volume curve and may be the signature of a liquid-liquid phase 
transition, in analogy with water |120j . Van Hoang has carried out a similar study 
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Figure 10. Upper panel: Simulated (dashed) Neutron structure factor T{r) 
compared to experimental findings |57j . Lower panels: Partial structure factors of 
300 K amorphous GeOi using the Oeffher-Elliot (solid line) and the Morse-like 

potential (broken line) |114j . 
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|121| under the same kind of conditions, i.e. in the high temperature (5000 K) liquid 
with changing density, and has suggested that a diffusion maximum should be attained 
for a density of about 5 g/cm 3 . This anomaly appears to be produced by competition 
between the breakdown of the tetrahedral network structure leading to an increase in 
atomic mobility, and the packing effects arising from densification that tend to reduce 
the mobility. The simulated structure of liquid Ge0 2 and Si0 2 appears to be very 
similar when the partial atomic correlation functions are properly rescaled [122 . 

Micoulaut, Guissani and Guillot [59] have used the Oeffner-Elliot potential to study 
the glass and liquid phases and which allows comparison with experiments. In the glass, 
the structural properties can be simulated relatively well, even though some structural 
limitations of the potential appear. While the first structural peak due to Ge — O 
interactions can be modelled very well at the expected distance of 1.72 A, as can the 
— distance at 2.81 A, the Ge — Ge correlations appear to be slightly overestimated 
(3.32 A) with respect to experimental values. This overestimation leads to a larger 
calculated value for the intertetrahedral angle than that obtained experimentally; 159° 
versus 130°, respectively. It is now well known that simple ionic potentials such as the 
ones reported above |l()9j . [TT] results in Ge — O — Ge angles that are too wide, a situation 
that has been encountered and reported already for amorphous silica [123 . However, 
the absence of any Ge — Ge interaction in the effective Oeffner-Elliot potential, except 
in the Coulombic term, may be responsible for the increased distortion in germania 
with respect to silica. In spite of these deficiencies, the simulation correctly describes 
the structure factor S(Q) and the partial structure factors Sij(Q) (Figure 5) and allows 
one to infer the origin of the first sharp diffraction peak (FSDP) as mostly arising from 
Ge — Ge correlations. However, overall the potential is found to reproduce the features 
of neutron scattering functions (Figure 10) reported by different groups [33], [5T] . 

Simulation using a Morse-like potential |114j provides a somewhat better agreement 
with the experimental partial atomic correlation functions (Figure 10) as the Ge — Ge 
distance is found to be 3.21 A at 300 K with correct bond angles [Qo-Ge-o — 108° and 
Oce-o-Ge = 133°) whereas both Ge — O and — distances are slightly underestimated 
(1.69 A and 2.78 A respectively) relative to experimentally derived values. 

6.3. Glass transition problem of strong glasses 

Enthalpy and glass transition temperature can be simulated rather well [SS] with 
respect to calorimetric measurements |124j . With the Oeffner-Elliot potential, a T g 
of 900 K, is found from the inflexion point of the potential energy. This value is close 
to the experimental derived T g (850 K, |125j ). This appears to be rather unusual 
as MD simulations on similar systems |126j . |127j predict much higher glass transition 
temperatures than the corresponding experimental ones. This is partially due to the high 
quench rates applied. In the present simulated systems, onset of slow dynamics at the 
nanosecond scale occurs in the same range of temperatures (920 K) which corroborate 
the calculated T g from the inflexion point of the energy profile. 
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Figure 11. Simulated diffusion for germanium and oxygen using the Oeffner-Elliot 
potential 

When put in contrast with silica, a more careful inspection of the self-diffusion 
coefficient D with respect to the viscosity behaviour |128j shows that the agreement 
between the simulated and experimentally measured T g reveals an underlying failure of 
the simulation technique. The self-diffusion coefficient D is computed from the mean- 
squared displacement of the germanium and oxygen atoms, and shows Arrhenius like 
behaviour D = D exp[£7j/T] at low temperatures, whereas at higher temperatures 
(T > 1600 K) some curvature appears (Figure 11), similar to that found for molten 
silica [129J. However, the calculated oxygen diffusion constant for GeC>2 at 1440 K 
is several orders of magnitude larger than the reported data for oxygen diffusion 
{Do = 7 x 10~ 14 m2.s~ 1 , [130J). A predicted diffusion D constant from viscosity 
data r\ using the Eyring relation k^T '/rjD = A (where A is a hopping length of about 
several Angstroms, |131| ) shows that both silica and germania overestimate the diffusion 
constants with respect to their simulated T 5 's thus allowing the system to remain in a 
liquid-like behaviour to lower temperatures. This underscores both the limitation of 
the employed potentials and the size of the simulated systems (actually up to several 
thousands atoms) to accurately describe the glass transition of strong glass formers. 

6.4- Equation of state 

The equation of state (EOS) of GeO<i has been reported by different authors, either 
on the basis of simulations |119j or from empirical models based on simple structural 
arguments [132|- In the latter, Smith and co-workers have shown that a two state 
function, taking into account the effect of the tetrahedral and octahedral character at 
low and high pressures, is able to describe the experimental equation of state at 300 K, 
whereas molecular dynamics simulations only succeed in simulating the EOS in the low 
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pressure range. 

At higher temperatures and higher densities, Gutierrez and Rogan |119j have shown 
that for simulated GeC>2 in the 3.5 — 5.6 g/cm 3 and T = 1500 — 3000 K range, 
pressure displays a monotonic decrease with molar volume. In the same context, a 
Birch-Murnhagan type |133j of EOS has been used [59J to fit a set of 269 simulated 
state points in the thermodynamic diagram. The method allows the extraction of the 
isothermal compressibility k,t as a function of temperature and density for density ranges 
lying between the ordinary glass density at 300 K (p = 3.66 g/cm 3 , jHIj) and about 
2.5 g/cm 3 . Progressive deviation of the Birch-Murnhagan EOS with respect to the 
simulated thermodynamic points appear for p < 2.5 g/cm 3 at high temperatures. At 
2000 K, the computed compressibility (n? = 9.13 x 10~ n Pa^ 1 ) is rather close to the 
experimentally measured value of Dingwell et al. (k,t = 12.4 x 10 -11 Pa -1 , |134j ). 

Micoulaut and Guissani (HH] have used a Direct Molecular Dynamics Method |135j 
to follow the equation of state at zero pressure, in order to predict the liquid-vapor 
coexistence curve of germania on the low (vapour) and high density (liquid) side in order 
to compare it with experimental results in the liquid up to 1440 K jflTj. Furthermore, the 
method highlights the quality of the effective potentials employed at low temperature. 
At zero pressure and low temperature (300 K) , the density of a simulated Oeffner-Elliot 
Ge0 2 glass |109j is indeed 3.70 g/cm 3 whereas the density of a simulated Ge0 2 glass 
using an alternative potential jTTj substantially disagrees with the experimental low 
temperature density of the liquid (p = 4.25 g/cm 3 as compared to the experimental 
p = 3.66 g/cm 3 ). Note however that this potential was used to study pressure induced 
rigidity in Ge0 2 (see below, |136p and the density at zero pressure was found |137j to 
be 3.9 g/cm 3 , i.e. much closer to the experimental value. The thermal history of the 
simulation appears therefore to be crucial in this case. 

Using a Wegner type expansion |138j , a critical point for germania is predicted and 
is located at T c = 3658 K, p c = 0.59 g/cm 3 and P c = 40 MP a [HI]. For the Tsuchiya 
potential jTTj, the location of the critical point seems to be much higher in temperature 
|14()j . This shift may arise from the increased charges used in the effective potential. 

6.5. Pressurized germania 

The application of pressure to amorphous germania seems to affect the structure 
stepwise. Experimentally a jump in bond distance from 1.72 A to 1.86 A is observed 
at around 9 GPa, signalling the conversion of tetrahedral to octahedral local structure 
as already described. However, numerical simulations show jBHj, at least in the low 
pressure range, that this conversion is somewhat more subtle. For pressures up to 
2 GPa, long-range correlations are reduced, as seen from the shift to higher wave vector 
of the position of the FSDP; similar to experimental observations |141j . In addition, 
a reduction is observed in the intertetrahedral bond angle (Ge — O — Ge) and then 
for P = 3 GPa a sharp distortion of the GeO^ tetrahedron occurs (Figure 12). These 
results are accompanied by a global increase in the number of oxygen neighbours in the 
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Figure 12. Local structure of germania with applied pressure . Right panel: 
Distortion parameter Sx of regular GeO^ and SiO± tetrahedra, as a function of applied 
pressure (X — Ge, open triangles). For comparison, the same parameter for S1O2 
(X = Si, filled triangles) is shown. The broken horizontal line represents the value of 
the perfect tetrahedron 5 = y/3/8. Right panel: Mean bond angles Ge — O — Ge and 
O — Ge — O with respect to compression (open and filled symbols) and decompression 
(broken curves). 

vicinity of a germanium atom that parallel the increase in density [59*] . 

Sharma and co-workers |142j have studied both the pressure induced structural 
changes of the a-quartz-like GeC>2 polymorph and amorphous GeC>2 using the Oeffner- 
Elliott potential in the (N,P,T) ensemble. The results show that both the average bond 
distance {Ge — O) and the average Ge coordination in a-quartz-like GeOi undergo a 
sharp change at around 8 GPa under compression, similar to the experimental findings 
of Itie et al. (H]. On decompression, the denser phase transforms back to a lower- 
density phase at ~ 2 GPa. The details of the number of oxygen neighbours around a 
Ge atom shows however, that the high density phase is not fully six-fold coordinated, 
as about 15% five- fold and 20% four-fold germanium can be found (Figure 13). Less 
abrupt changes are expected for vitreous GeO^ (Figure 13) where a majority of six-fold 
germanium only occurs for pressures larger than 20 GPa. 

The structural changes with pressure are more dramatic in liquid (1650 K) GeOi as 
a sudden loss of five-fold germanium atoms and an almost six-fold coordinated structure 
is obtained for pressures larger than 12 GPa. 

Finally, it appears that the evolution of the intermediate range order with pressure 
or density is selective as MD simulated ring statistics [143^ show that rings with more 
than six germania tetrahdera tend to disappear for densities larger than 5 g/cm? whereas 
the growth of edge-sharing GeO^ octahedra signals a behaviour similar to Ti02- 
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Figure 13. Calculated variation with pressure |142| of the fractional abundance of Ge 
coordination in the a-quartz-like GeOi polymorph, glassy (300 K) and liquid GeC>2 
(1650 K). 

6. 6. Pressure induced rigidity and intermediate phases 

Trachenko et al. |144j . ^45] have been investigating the network rigidity of GeC>2 
and S1O2 under pressure. Rigidity usually appears when the number of mechanical 
constraints per atom, arising from interatomic interaction (mostly bond stretching and 
bond bending) becomes greater than the number of degrees of freedom |140j . In network 
glasses, this is generally achieved by the addition of cross-linking elements such as 
germanium into a basic flexible structure containing e.g. selenium chains. This leads 
to an increase of the network mean coordination number f (and to the increase of 
constraints) and produces a stiffening of the structure and ultimately a floppy to rigid 
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Figure 14. Relative variation in volume in GeO-i and S1O2 glasses under pressure 
compared with experiments on SiC>2 

transition. The onset of rigidity and the way it percolates has been documented for 
various glass- forming systems. In recent years however, a reversibility window |147j . 
|148j has been discovered located between the floppy and rigid phases, and which 
manifests itself by the loss of irreversibility (and hysteresis) of the heat flow when cycling 
through the glass transition temperature region. A similar state can be found in glassy 
GeC>2 and SiC>2 under pressure. 

Pressure induced rigidity in Ge0 2 glass using MD simulations has been addressed 
recently |13fij . Here the increase in connectivity (or mean coordination number f) 
is achieved with the increase of the glass density or the application of pressure that 
produces a tetrahedral to octahedral conversion. Thus pressure introduces locally 
rigid higher-coordinated units in an otherwise flexible tetrahedral network of GeO^ 
tetrahedra. Densification with temperature under pressure can take place in a pressure 
window, centred around the rigidity percolation transition. The density change is about 
7%. This new effect has been rather well documented for silica |144j . |145j and compared 
successfully with experimental results, |149j and additional simulations showing the loss 
of low frequency modes in the effective vibrational densities of states at the same pressure 
where densification occurs. However, it has only been shown that the pressure window 
in germania is centred around 0.5 GPa, i.e. considerably lower than for silica (5 GPa) 
(Figure 14). However, it is another signature of the increased sensitivity to pressure 
change of Ge02 with respect to Si02- Indeed, the tetrahedral to octahedral conversion 
of amorphous S1O2 manifested by the jump in Si — O bond distances is found to be 
around 13 GPa |15Uj whereas the same jump is found to be at 8 GPa for Ge0 2 glass 
|Hj. It is therefore not surprising at all that onset of rigidity manifests at lower pressures 
in Ge0 2 - 
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6. 7. Ab initio studies of c-GeC>2 and germania 

One way to circumvent the possible failures of the above mentioned semi-empirical 
potentials, is the use of ab initio methods, especially under extreme conditions where 
the potentials are not necessarily reliable. Hafner and co-workers have studied the 
high-pressure transformations up to 70 GPa of crystalline GeC>2, using density functional 
theory with a pseudopotential method, and a local density approximation. It appears 
from this computation that several high-pressure phases can exist in GeC>2 which are 
a tetragonal CaCl 2 type at 40 GPa, an alpha — Pb0 2 -type at 40 GPa and finally 
a pyrite-type crystal at 70 GPa, similar to those observed experimentally (see above) 
(Figure 3). These transformations highlight the analogy of the phase transition sequence 
between Si0 2 and Ge0 2 polymorphs at high pressure. Additional studies concerning the 
electronic properties of these polymorphs have been reported by Christie et al. |151j . 
using the same tools. This allows determination of the lattice parameters, cohesive 
energy and bulk modulus by minimizing the total energy of the solid. In addition, an 
equation of state for the polymorphs can be fitted with a Birch-Murnhagan EOS |133j . 
|152j or the density of states. 

Ab initio studies of amorphous germania have been only reported recently [SB! using 
the same numerical scheme but with an improved density approximation (generalized 
gradient). This enables determination of the neutron structure factor (Figure 4), the 
infrared and the Raman spectra, all of which show good agreement with experimentally 
derived data (cf., Figure 8c). It furthermore provides insight into the 3-membered 
ring distribution and the so-called D 2 line first described by Galeener and workers [82J. 
The projection of vibrational eigenmodes onto natural or isotopic substituted oxygen 
breathing motions in these rings shows that a broad peak centered at 520 cm -1 and 
corresponding to the experiment, is blueshifted with 18 by 26 cm -1 . The number of 
these rings is found to be about 20% of the oxygen atoms. On the other hand, similar 
calculations do not seem to support the assignment of four-membered rings to the Di 
line found at 347 cm' 1 . Instead, this band seems to arise from coupled motions of Ge 
and O atoms. 

7. Summary and Conclusions 

Studies on the structure of crystalline, liquid and glassy Ge0 2 continue to be of interest 
to a number of researchers in physics and glass, materials and geological sciences. This 
breadth of interest stems from the fact that while there are close similarities between 
Ge0 2 and Si0 2 , there are also distinct differences which make Ge0 2 useful as an 
analogue for studying the high pressure behaviour of oxide glasses. Crystalline Ge0 2 
polymorphs behave, with increasing temperature and pressure, in a manner similar 
to crystalline Si0 2 polymorphs. However, pressure induced phase transformations 
generally occur at much lower pressures than equivalent Si0 2 phases. This is because 
the larger Ge0 4 tetrahedron (relative to S1O4 tetrahedron) is more distorted due to 
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greater variability in the O — Ge — O angles. This makes the use of Ge0 2 polymorphs 
attractive as Si0 2 analogues in high pressure studies for studying possible pressure 
induced structural changes, since the pressure ranges required are much more accessible. 

Ge0 2 glass has also been considered as being somewhat similar to Si0 2 glass. The 
first three interatomic distances in the glass are reasonably well resolved and indicate 
that, like Si0 2 glass, the network is composed of tetrahedra linked together through 
their corner bridging oxygens. However, there are significant differences between the 
two glass networks. Ge0 2 glass has a much smaller mean Ge — O — Ge angle and 
a much higher proportion of 3-membered rings, relative to Si0 2 glass. Furthermore, 
there may be differences in the intermediate-range structure with Ge0 2 glass possibly 
being composed of 4- membered rings, rather than the currently accepted 6-membered 
rings similar to Si0 2 glass. In addition, application of high pressure readily converts 
4-fold Ge to 6-fold Ge, via a transitionary 5-fold coordination, at much lower pressures 
than found for Si0 2 glass. In the liquid state Ge0 2 retains 4-fold geometry to high 
temperature but with broadened Ge — O — Ge angles although numerical studies of 
liquid Ge02 indicate that while the Ge0 4 tetrahedra are distorted, the Ge — O — Ge 
angle remains similar to that found in the a-quartz-like Ge0 2 polymorph. Furthermore, 
simulation of pressure effects indicate that the pressure-induced transformation from 4- 
to 6-fold Ge observed experimentally at ~ 9 GPa, may be quite subtle. Onset of the 
simulated transition occurs at 2 GPa with loss of long-range correlations, a reduction in 
the Ge — O — Ge angle, followed by a sharp distortion of the GeO^ tetrahedra at 3 GPa. 
This is accompanied by the onset of rigidity at much lower pressure than observed for 
Si0 2 glass. When Ge substitutes for Si along the Ge0 2 — Si0 2 binary, there is no 
evidence for clustering or phase separation of the glass network and it is composed of 
SiO/i and GeO^ tetrahedra. The substitution is random with no heterogeneity induced in 
the combined network. However, with increasing pressure Ge undergoes a coordination 
change from 4- to 6-fold coordination. The pressure at which this occurs is dependent 
upon the Si0 2 composition indicating that Si has an influence on the local structure of 
Ge. In addition, there is a broad pressure-composition range over which Ge is in both 
4- and 6-fold coordination. 
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